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ABSTRACT 

High-redshift, dust-obscured galaxies - selected to be luminous in the radio but relatively faint at 
850 /im - appear to represent a different population from the ultra-luminous submillimeter- (submm-) 
bright population. They may be star-forming galaxies with hotter dust temperatures or they may have 
lower far-infrared luminosities and larger contributions from obscured active galactic nuclei (AGN). 
Here we present observations of three z ~ 2 examples of this population, which we term submm- 
faint radio galaxies - SFRGs in CO(3-2) using the IRAM Plateau de Bure Interferometer to study 
their gas and dynamical properties. We estimate the molecular gas mass in each of the three SFRGs 
(8.3 X IO^Mq, < 5.6 X IO^Mq and 15.4 x IO^Mq, respectively) and, in the case of RG J163655, 
a dynamical mass by measurement of the width of the CO(3-2) line (8 x 10^" csc^ zM©). While 
these gas masses are substantial, on average they are 4 x lower than submm-selected galaxies (SMGs) . 
Radio-inferred star formation rates (< SFRiadio >= 970 Mq yr~^) suggest much higher star-formation 
efficiencies than are found for SMGs, and shorter gas depletion time scales (~llMyr), much shorter 
than the time required to form their current stellar masses (~160Myr; ^10^^ M0). By contrast, SFRs 
may be overestimated by factors of a few, bringing the efficiencies in line with those typically measured 
for other ultraluminous star-forming galaxies and suggesting SFRGs are more like ultraviolet- (UV- 
)selected star-forming galaxies with enhanced radio emission. A tentative detection of RG J163655 
at 350 /im suggests hotter dust temperatures - and thus similar gas-to-dust mass fractions - as the 
SMGs. We conclude that SFRGs' radio luminosities are larger than would naturally scale from local 
ULIRGs given their gas masses or gas fractions. 

Subject headings: galaxies: evolution — galaxies: formation — cosmology: observations — galaxies: 
starbursts — galaxies: high-redshift 



1. INTRODUCTION 

Submm surveys have provided an efficient probe of 
star-formation activity in ultraluminous infrared (IR) 
galaxies (ULIRGs, >1O^^L0) in the distant Universe 
(e.g., Small, Ivison & Blain 1997; Hughes et al. 1998; 
Barger et al. 1998), with bright submm emission pro- 
viding unambiguous evidence of massive quantities of 
dust, heated predominantly by young stars rather than 
AGN (e.g.. Chapman et al. 2003a; Alexander et al. 2005; 
Menendez-Delmestre et al. 2007; Valiante et al. 2007; 
Pope et al. 2008). Before the availability of the Ata- 
cama Large Millimeter Array (ALMA), confusion will 
continue to limit the sensitivity of current submm sur- 
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veys. As a result, many ULIRGs fall below the detection 
limits due to variations in their spectral energy distribu- 
tions (SEDs) - usually parameterized in terms of dust 
temperature (Td) - meaning that entire populations of 
star-forming galaxies may have been missed by submm 
surveys. 

For a fixed far- infrared luminosity (FIR) , a galaxy with 
a higher Td will be weaker in the submm at 850/im than a 
galaxy with a lower T^- Specifically, raising Td from the 
canonical ^^35 K for SMGs to 45 K will result in a factor 
'--^lOx drop in 850/im flux density (Blain 1999; Chapman 
et al. 2004). These galaxies should, though, be accessible 
in the radio waveband, regardless of their specific SEDs, 
since the radio correlates with the integrated FIR emis- 
sion (Helou et al. 1985) with a small ~ 0.2dex dispersion 
and no observable dependence on SED type. However, 
there is potential for large AGN contaminations in the 
radio, as has often been the case with mid-IR selection of 
z> 1 ULIRGs (e.g., Houck et al. 2005; Yan et al. 2005, 
2007; Sajina et al. 2007; Weedman et al. 2006a, 2006b; 
Desai et al. 2006) and the facilities required to provide 
the high-resolution, multi-frequency radio data needed 
to decontaminate the samples (e.g., Ivison et al. 2007a) 
are not yet available. 

Substantial populations of apparently star-forming 
galaxies at z ~ 2 have been uncovered through deep 
1.4-GHz radio continuum observations, many of which 
are not detected at submm wavelengths with the current 
generation of instruments (Barger et al. 2000; Chapman 
et al. 2001, 2003b, 2004a - hereafter C04). These galax- 
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ies are luminous in the radio and spectroscopy suggests 
that star formation is powering their bolometric out- 
put (there is httle or no sign of high-ionization emission 
Hues, characteristic of AGN in their UV/optical spec- 
tra). These galaxies could, in principle, span a range in 
properties from deeply obscured AGN to far-IR-luminous 
starbursts. In the latter case, one would expect a differ- 
ent SED from a typical SMG - a higher T^ for instance. 
These submm-faint radio sources have a large volume 
density at z '^ 2, even larger than the SMGs (Haarsma 
et al. 1998; Richards et al. 1999; Chapman et al. 2003a, 
C04; Barger et al. 2007). There are p == 2 x 10"^ Mpc^ " 



radio sources with Li.4ghz > 



10^^ ergss ^ 



Hz- 



at z 



compared with p = (6.2 ± 2.3) x 10"^ Mpc"^ foj. sMGs 
brighter than 5 mJy at 850 /im at the same epoch (Chap- 
man et al. 2003b, 2005). As essentially aU of these SMGs 
form a subset of these radio sources (Chapman et al. 
2005; Pope et al. 2006; cf. Ivison et al. 2002) this im- 
plies '^ 14 X 10^^ Mpc^"^ luminous radio sources remain 
undetected at submm wavelengths. Understanding the 
exact properties of these galaxies is therefore of great 
importance. If they are all forming stars at the rates im- 
plied by their radio luminosities, they would triple the 
observed SFR density (SFRD) at z ~ 2. By contrast, if 
their radio luminosity comes from a mix of star formation 
and AGN, they have less impact on the global SFRD but 
they increase the highly obscured AGN fraction at these 
epochs (e.g., Daddi et al. 2007a; Casey et al. 2008) and 
contribute substantially to black hole growth. 

Together with other observations, the redshifted cool- 
ing emission lines of CO allow us to assess and compare 
the energy source of SFRGs with that of SMGs and other 
distant star-forming galaxies via measurements of their 
gas and dynamical masses. In this paper, we present 
the results of a pilot study with the IRAM - Plateau de 
Bure Interferometer to detect molecular gas in SFRGs 
through the rotational CO (3-2) line emission. In § 2 
we describe the sample properties and observations both 
with PdBI and other facilities. Section § 3 presents 
the CO(3-2) detections and limits obtained from the 
PdBI observations, § 4.1 estimates gas properties, star- 
formation rates and efficiencies, and § 4.2 compares the 
SFRGs to other galaxy populations. Finally § 5 discusses 
the results and places them in a broader galaxy evolu- 
tion context. Throughout we assume a cosmology with 
h = 0.7, f^A = 0.72, r^M = 0.28 (e.g., Hinshaw et al. 
2008). 

2. SAMPLE PROPERTIES AND OBSERVATIONS 

Our sample is drawn from an expansion of the C04 
submm-faint, radio-selected galaxy (SFRG) program, 
with galaxies drawn from several deep radio survey fields 
with typical sensitivity limits of ct = 4 — 8yLtJy (e.g., Biggs 
& Ivison 2006). In the submm, the survey fields are im- 
aged to a typical depth of crsso/jm ~ 1 — 2mJy (e.g., 
Scott et al. 2002, Borys et al. 2003). We selected sources 
with redshifts, radio luminosities, and submm limits typ- 
ical of the population (< z >= 2.1, < Liaghz >= 
2 X 10^1 ergss-^Hz-i, Lgso^m < 1 x lO^^ ergss-^Hz"! 
(< 2 mJy for z ~ 2) at the ~ 2ct level, lying within ±1ct 
of the median SFRG in (C04), and observing RG J163655 
and RGJ131236 based on suitability of RA and con- 
fidence in the optical spectroscopic redshifts. A third 
source from our sample, RG J123711, was observed pre- 
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Fig. 1. — top panels: CO(3-2) spectra for the two candidate de- 
tections. The spectra are shown smoothed with a 50 km/s boxcar 
filter, and with respect to the zero velocity offsets defined from the 
ffa emission line redshift (red dashed line). The best- fit Gaussian 
profile is shown for the emission line in RG Jf 63655, along with the 
Uv-fnferred redshift from inter-stellar absorption lines (blue dashed 
line), bottom panels: velocity-averaged spatial maps of CO emis- 
sion, from —1500 to — 800kms~^ (RG J163655) where contours are 
from —1 to 5(7 in steps of a (0.05mjy beam~^). RG J123711 does 
not represent a formal CO detection. We measure a significance of 
3.2(T integrating over the full band). The field of view, on a side, is 
35", with the size of the beam shown to the lower left. Both CO 
emitters lie exactly at the radio source position to within 1". 
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viously with PdBI in the same field as an SMG in the 
program described by Neri et al. (2003) and Greve et al. 
(2005), and we include this object here. We note that 
while the neighboring SMG (SMM J123712.0+621326) 
lies only 8" to the south-east, we are confident that 
RG J123711 is not a luminous submm emitter. Firstly, 
SMM J123712, has a strong CO line detection (Small et 
al., in preparation) with a ~ Scr detection 5*00(3-2) = 
1.2Jykms~^, comparable to typical SMGs from Greve 
et al. (2005). Secondly, the 850-/^m emission peak (with 
a R=7" beam) in our SCUBA map is centered on the 
radio source position, whereas no significant peak is ob- 
served at the position of RG J123711. Removing a 850- 
/im point source from the position of SMMJ123712 re- 
veals an even lower 850-/im flux density (0.1 ± 1.2 mJy) 
at the position of RG J123711 than the 2.2±1.2mJy con- 
servatively adopted for our calculations (Table 2). Im- 
portantly, the small 850-/im/1.4-GHz flux ratio for this 
source is clearly comparable to SFRGs and not to the 
typical radio-detected SMGs in Chapman et al. (2005) 
or Ivison et al. (2007b). The properties of these SFRGs 
are listed in Tables 1 and 2, and displayed in Figs 1 and 
2. 

2.1. PdBI observations 

RGJ163655 and RGJ131236 were observed in their 
redshifted CO (3-2) lines and in the continuum at 
'^lOSGHz using the newly refurbished PdBI receivers 
for 11.5 and 4.8 hr, respectively. Observations were 
made in D configuration on 2007 January 24, April 
28, May 08 and June 03, with good atmospheric phase 
stability (seeing, 0.7-1.4") and reasonable transparency 
(0.5mm of precipitable water vapor). For RG J163655, 
the Ha redshift showed a considerable offset from that 
inferred from UV absorption and emission lines. While 
this was conceivably due to a large velocity starburst 
wind (~1300kms~^), we considered the possibility that 
one redshift might have calibration or resolution prob- 
lems. We therefore observed RG J163655 split over two 
slightly offset frequency settings to span all measured 
redshifts. The overall flux scale for each observing epoch 
was calibrated using a variety of sources. In each ob- 
serving epoch between three and six sources were used. 
The visibilities were resampled to a velocity resolution 
of 55kms~^ (20 MHz) providing la line sensitivities of 
1.6mJy beam~^. The corresponding synthesized beam, 
adopting natural weighting, was similar for both sources, 
5.0" by 4.0" at PA '-SO degrees, east of north. Obser- 
vations of RG J123711 proceeded similarly to those de- 
scribed by Greve et al. (2005). The PdBI data for aU 
three SFRGs were calibrated, mapped and analyzed us- 
ing the GILDAS software package. The CO (3-2) spec- 
tra and images of the two candidate SFRG detections, 
RG J163655 and RG J123711, are shown in Fig. 1. 

2.2. Other Observations 

Spitzer observations of the SFRGs in this paper were 
taken with the Infrared Array Camera (IRAC) and the 
Multi-band Imaging Photometer (MIPS) through various 
GTO and Legacy programs. 

RG J163655: In addition to the SCUBA 850-Aim pho- 
tometry, observations were obtained at 350 /im with the 
SHARC2 camera (Dowell et al. 2005) on the Caltech Sub- 
millimeter Observatory as part of the imaging campaign 
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Fig. 2.— SEDs for RG J163655, RG J131236 and RG J123711. A 
local galaxy SED with Tj = 40 K (from Dale et al. 2003) is shown, 
normalized to the radio point in each case. An average fit to the 
SEDs of representative SMGs is shown as a heavy line (Chapman 
et al. 2005; Kovacs et al. 2006; Hainline et al. 2008). 



of Kovacs et al. (2006), where observational details can be 
found. A near-IR spectrum from UKIRT/UIST covering 
the Ha/[Nii] region (Swinbank et al. 2006) finds a large 
Ha/ [N ii] ratio suggestive of a relatively low metallicity, 
[Fe/H] r^ —0.9, and no strong AGN component. While 
the [Oiii]A500.7 region was also covered with the instru- 
ment, we did not clearly detect any emission lines, setting 
a limit on the [Oiii]A500.7 line of 3.75 x 10"" Wm'^ 
again suggesting that an AGN does not dominate the 
energetics of this galaxy. The line width of the Ha emis- 
sion (FWHMrost) is 420±100kms-i and the integrated 
line flux of 5.0 x 10~^^ W m~^ then suggests a SFR uncor- 
rected for dust extinction of 15O±5OM0 yr~^ (Kennicutt 
1998). The 1.4-GHz radio emission was unresolved by the 
Very Large Array (VLA) in its A configuration (Biggs & 
Ivison 2006). A relatively compact UV morphology is 
observed in the Hubble Space Telescope imaging, with 
Ri/2=0.25" (Swinbank et al. 2006). 

RG J131236: In C04, the optical (rest-frame UV) spec- 
trum was presented, showing Lya in emission but all 
other detectable lines in absorption, and was classified as 
a pure starburst. The near-IR spectrum of RG J131236 
from Keck/NIRSPEC covering the Ha/[Nii] region 
(Swinbank et al. 2004), again finds a large Ha/ [Nil] 
ratio suggestive of a low metallicity, [Fe/H] ~ —0.9. 
The line width of the Ha emission is 450 ± 220kms~^ 
and the integrated line flux of 2.3 ± 1.2 x 10~^^ Wm^^ 
then suggests an SFR uncorrected for dust extinction of 
110 ± 40MQyr~-'^. The radio emission was unresolved 
by the VLA. Only ground-based imaging exists for this 
SFRG, showing a faint, unresolved source in 0.8" seeing. 

RGJ123711: An X-ray-detection and obscured AGN 
classification due to its X-ray luminosity by Alexander et 
al. (2005) contrasts a PAH-dominated mid-IR spectrum 
obtained by Pope et al. (2008), showing no obvious AGN 
component. No rest-frame UV spectrum for RG J123711 
has been published as there are no detected features. The 
redshift is based entirely on the near-IR spectrum where 
RGJ123711 was detected with Keck/NIRSPEC (Swin- 
bank et al. 2004) , once again finding a large Ha/ [N ii] 
ratio suggestive of an [Fe/H] ~ —0.8. The narrow line 
width of the Ha emission (112 ± 45kms~"'^) and the in- 
tegrated line flux of 0.4 ±0.3 x 10~^®Wm~^ suggests an 
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SFR uncorrected for dust extinctfon of 16 ± 9 Mq yr ^. 
Radio imaging of this SFRG with MERLIN (0.3" synthe- 
sized beam; Casey et al. 2008b) reveals a double source 
structure with a ~1" elongated feature and a relatively 
compact Ri/2=0.4" component. 

3. RESULTS 

The velocity-integrated line fluxes or limits for all three 
SFRGs are hsted in Table 1. For RG J163655, inspection 
of the data cube shows a significant 4.9cr detection of 
CO (3-2) line emission at the phase center, integrated 
over the velocity channels at ^ — lOOOkms"^ with a ve- 
locity width of 400kms~^ fwhm. Fitting a Gaussian 
profile to the CO spectrum, we derive a best-fit redshift 
for the CO(3-2) emission of z = 2.1859 ± 0.0002, and es- 
timate the CO fiux is by summing the channels from —2a 
to +2(7 of Gaussian fit to the line. We note that no sig- 
nificant continuum emission is detected from the line-free 
region (^-^650 MHz of bandwidth) down to a la sensitiv- 
ity of 0.07 mjy beam^^ , consistent with the submm limit, 
assuming a dust spectral index, v'^^-^, for a modified 
blackbody with emissivity, (3 = +1.5. For RGJ131236 
no significant emission is observed at the phase center, 
although the limit on the CO gas mass is still of great 
interest relative to the SMGs. For RGJ123711 we ten- 
tatively detected (3.2cr) a positive signal integrated over 
the full band from —600 to +600 km s^^, centered on the 
Ha determined redshift of 2; = 1.996. A precise CO red- 
shift cannot be determined for RGJ123711 as the line 
shape cannot be determined in the low S/N spectrum. 
To assess the possibility that we have simply detected 
continuum in this source, we analyze the radio spectral 
index, measured to be steep from the 8.4 GHz/1.4 GHz 



flux ratio (Muxlow et al. 2005), S,. 



, and the 



synchrotron contribution at ^^3 mm is negligible. This 
does not however preclude a contribution from an AGN 
component with the opposite spectral slope emerging at 
higher frequencies. It is clear that the full radio spectrum 
is needed to explore this issue, as well as the possibility 
of an obscured AGN. 

In RG J163655, the CO-inferred redshift is close to the 
redshift measured from various interstellar absorption 
lines in the Keck/LRIS UV spectrum, but is blue shifted 
by llOOkms^^ from the Ha line detected by Swinbank 
et al. (2006). Re-analysis of the near-IR spectrum does 
not signiflcantly change the result as the sky line cali- 
brations appear to be as presented in Swinbank ct al. 
(2006). If the detected line were dominated by [Nil] with 
Ha/[Nii]< 1 - and we stress that there is no evidence 
in the spectrum of this - then the implied velocity off- 
set would be ^^800 km s^^, somewhat closer to the aver- 
age CO velocity and consistent with the higher velocity 
peak in the detected CO profile. We cannot attribute the 
CO emission to an offset companion (as was the case for 
SMC J09431 in Tacconi et al. 2006) as the CO centroid 
is exactly at the near-IR and radio position to within 
the 1" centroiding uncertainty (^^beam size x (S/N)^^). 
The Ha-inferred redshift has not always been represen- 
tative of the CO redshift in SMG surveys (the average 
CO-Ha offset is 150 kms"-'^), presumably because the lu- 
minous core is so deeply dust enshrouded that wind out- 
flows or satellite H 11 regions are more strongly detected 
in Ha. We therefore put forward the hypothesis that we 
are detecting a highly dust-obscured gas-rich galaxy in 



CO (3-2), either a companion seen in projection or else 
one not well sampled by the Ha observations. The rest- 
frame FWHM of the CO line is 376 ± 40kms"^, close to 
that found in Ha by Swinbank et al. (2006), but likely 
a coincidence given that the CO and Ha redshifts are 
discrepant. 

The 350Aim SHARC-2 imaging of RG J 163655 shows a 
tentative continuum detection, Saso^jm — 2.4±6.5mJy at 
the radio position, however given the telescope pointing 
errors and the low signal-to-noise (S/N) of any expected 
emission, we search a region comparable to the beam size 
(9"). A 2a peak (S'sso^.m = 13.7 ± 6.9 mJy) lies 7" from 
the radio position at 16h 36m 54.6s, +41° 04' 28" J2000. 
Within this area there are only ^^3 SHARC-2 beams and 
the chance of a spurious 2a peak is only ~6%. There 
is a high likelihood (^90%) that this peak is related to 
RG J163655 and this flux range is completely consistent 
with our SED flt to the radio photometry for RG J 163655 
(Fig. 2). 

We calculate from the integrated CO (3-2) intensity 
(Jykms~^) the line luminosities and estimate the total 
cold gas masses (H2 + He) (listed in Table 1). We adopt 
^co(3-2) = 3.25 X 10^Sco{3-2)i^2i'^ + z)-^Dl, consis- 
tent with Solomon & Vanden Bout (2005). 

We assume both a line luminosity ratio of r^2 = 

co(3-2) _ ^ /J g^ ^ constant brightness tempcr- 

CO(l-O) 

ature) and a CO-to-H2 conversion factor of a = 
O.SM© Kkms"'^ pc^. These values are appropriate for 
local galaxy populations exhibiting similar levels of star- 
formation activity to our SFRGs (e.g., local ULIRGs - 
Solomon et al. 1997), and this choice also facilitates com- 
parison with SMGs modeled with the same values (Greve 
et al. 2005). We discuss later how the affect of adopting 
typical Milky Way values aco and r32. 

All three SFRGs show a peaked SED in the mid-IR 
(Fig. 2) suggesting that this spectral region is dominated 
by stars rather than AGN. These properties are used 
to derive the rest-frame K-hanA (~ 2.2/xm) fiux, and 
convert to a stellar mass, in a matter similar to Borys 
et al. (2005), adopting their Li<-/M0 = 3.2 characteristic 
of a burst with an age of ~ 250 Myr. We interpolated 
between IRAC bands to estimate 52.2^m (Table 3). 

3.1. Derived Properties 

We then proceed to estimate various derived properties 
(listed in Table 3). Starting with the gas surface den- 
sity, for RG J123711, we assume the CO emission traces 
the same large, extended morphology (> 1" FWHM di- 
ameter) traced by resolved MERLIN radio imaging (see 
Casey et al. 2008b for details), suggesting a low gas sur- 
face density. For RGJ163655 and RGJ131236, neither 
the CO emission nor the VLA radio emission is resolved. 
Without further information beyond the optical imag- 
ing described previously, we assume the gas in these two 
SFRGs is distributed in a disk with a similar radius to 
SMGs (e.g., Tacconi et al. 2008) of Ri/2 = 1.7kpc (0.25"), 
resulting in higher inferred gas surface densities The CO 
luminosities for the three SFRGs are plotted as a function 
of FIR luminosity and compared to other high-rcdshift 
galaxies detected in CO in Fig. 3. 

A dynamical mass for the well-detected RGJ163655 
can be estimated by analyzing the CO line profile. We 
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base our analysis on a single Gaussian fit to the line. CO 
emission is comparatively immune to the effects of obscu- 
ration and outflows and therefore provides a unbiased 
measurement of dynamics within the CO-emitting re- 
gion. The line width of the CO emission (410±40kms~^ 
implies a dynamical mass of (8.4 ± 2.1) x 10"'^° csc^ i Mq, 
assuming the gas lies in a disk with inclination i and a 
radius of 0.25" (1.7 kpc). Based on this, we calculate a 



gas to dynamical mass fraction of / 



M„ 



M. 



O.lOsin^i. 



We note that the mean angle of randomly oriented disks 
with respect to the sky plane in three dimensions is i 
~ 30° (Carilli & Wang 2006), resulting in an average 
inclination correction of csc^ z = 4. 

4. ANALYSIS 

4.1. Star- Formation Rate and Efficiency 

The radio luminosity of the SFRGs forms our base- 
line estimate for the far-IR luminosity and SFRs (listed 
in Table 3), since we have only upper limits at 450 and 
850 ^m. We caution however that our flux- limited ra- 
dio selection biases our sample to find objects of the 
same radio luminosity as SMGs, regardless of the ori- 
gin of the radio power. There are clear examples within 
the wider SFRG sample where AGN dominate the ra- 
dio power despite an apparent starburst spectrum in the 
rest- frame UV (e.g., Casey et al. 2008a). The average 

< SFRradio > = 97OM0yr"\ assuming the radio/FIR 
relation q = log(FIR/3.75 x 1012Hz)/S'i.4GHz) (Helou et 
al. 1985), with q = 2.34 (Yun et al. 2001), a correction 
factor of 2.3 to total IR luminosity (TIR) appropriate for 
hotter dust SEDs (Dale & Helou 2002), and the conver- 
sion from Kcnnicutt (1998) 

SFR(Af0yr-i) - 1.8 x lO-i"L8-ioooMm(i0). 

The SFRs from their dust-corrected rest-frame 1500A 
continuum flux are factors ^ 50 x less (as described 
in C04), and the UV emission is clearly not probing 
the true luminosities of these systems (Table 3). The 
Ha emission line suggests SFRs (Table 3) ten times 
less than the radio (< SFRhq >= 92), although with 
the average extinction factor of Ay ~ 2.9 ± 0.5 pro- 
posed for SMGs in Takata et al. (2006), this becomes 

< SFRhq^coit > = 1300 Moyr^^ This is consistent, 
on average, with the radio-inferred SFRs, although the 
individual radio/Hacorr ratios show very poor corre- 
spondence. We have, of course, argued in the case of 
RG J163655 that the CO(3-2) and the Ha line emission 
may be coming from distinct regions, so these arguments 
do not obviously apply in every case, and the average 
correction factor applied to the Lhq may not be appro- 
priate either individually or for the population. The 24- 
yum fluxes (Table 2) would represent strong supporting 
evidence for large SFRs. We calculate SFR24^m as in 
Pope et al. (2006) for consistency with SMGs (although 
strong 24-/im luminosities could also reveal a dominant 
hot AGN dust torus). In a large sample of SFRGs, the 
24-/im luminosity distribution is indistinguishable from 
that of SMGs (Casey et al., 2008b). Two of our present 
three SFRGs have 24-/im observations, and only one is 
detected, the limit in the second case not being partic- 
ularly constraining relative to the radio. The SFR24;^m 
(Table 3) could be consistent with the SFRradio given 
uncertainties in calibrating the SFR indicators. 



For a reference point - which can be scaled through 
by uncertainties in the SFR - we adopt the SFRradio 
and calculate surface densities Ssfr, and star- formation 
efBciencies (SFE — Lpm/MHa) for the SFRGs, using as- 
sumed sizes described previously, and listed in Table 3. 
We find large SFEs (4x larger than the SMGs on av- 
erage), although the observational constraints on the 
SFRs for the SFRGs are consistent with having been 
over-estimated by a factor ~ 2 — 3 x which would bring 
them into reasonable agreement with the envelope of 
star-formation efficiencies for SMGs, LBGs and local 
ULIRGs. 

We can also roughly estimate a lower limit to the gas 
to dust mass ratio where CO is detected (Table 3). We 
estimate similar dust mass limits for the SFRGs using 



M, 



dust 



"Mml = 



1 



SobsD 



1 + z KdB{vQ,Td) 



assuming k^ oc v^ , j3 = +1.5, and B^{Td) ^ v^'i, in the 
Mdust < 1-6 X 10^ Mq from the 2a limit on the 850 ^m 
flux of ^ 2.2 mjy, assuming a dust mass absorption co- 
efficient of K850 „ — 0.15 m^kg^ . We find ratios of 
50 and 100, with at least a factor of ~ 6 uncertainty 
accounting for our uncertainty of the dust temperature 
(ATd ~ ±5K), dust emissivity coefficient (A/3 ~ ±0.5), 
and mass absorption coefficient (about a factor of ~ 3; 
e.g. [Seaquist et al. 2004). 

Assuming the molecular gas reservoirs we detect are 
fueling the star formation within these galaxies, then 
there is enough gas to sustain the current star forma- 
tion for Tdopiction ~ M(H2)/SFR ranging from less than 
9 to 13Myr. Since we have also estimated stellar masses, 
we can compare the gas depletion time with the time 
to form the current stellar mass of the system. At the 
current SFRs, Tformation ~Mstars/SFR ranges from 80- 
240 Myr, which are comparable to the assumed ages of 
the stellar populations. 

4.2. Comparison to other populations 

Comparison of the SFRGs to the SMGs is of primary 
importance, since a major goal of the observations is 
to understand the degree to which SFRGs should be 
treated on a similar footing to SMGs in models and evo- 
lutionary calculations. Taking the average CO line lu- 
minosity and gas mass we find intrinsic line luminosi- 
ties a factor of ~4 lower than the median for SMGs 
(cf. < L'co >= (3.8 ± 2.3) X 10^° Lq a nd < A/gas >= 
(3.0 ± 1.6) X 10^° M^: lGreveeraII[200l . The CO line 
width of RG J163655 is also much lower than the median 
of SMGs (< FWHM >= 780 kms^^, Grevc et al. 2005). 
Given the CO-inferred gas masses are somewhat low 
compared to SMGs we find, not surprisingly, that the gas 
depletion timescales are short compared to SMGs (which 
have Tdcpiotion ~ 40 Myr, Greve et al. 2005), though the 
SFRGs are still within a physically plausible range. We 
only have a useful constraint on the size of the emitting 
region for RGJ123711 to compare with higher resolu- 
tion images of SMGs, although in general the SFRGs ex- 
hibit similar radio sizes and morphologies to SMGs (com- 
pare Chapman et al. 2004b and Biggs & Ivison 2008 to 
Casey et al. 2008b). The typical gas to dynamical mass 
fraction in SM Gs is estima.ted t o be ~ 0.3 assuming a 
merger model (jGreve et al.ll2005l) . while they have SFEs 
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Table 1 
Properties of CO observed SFRGs 



RA Dec z(CO) z(UV) z(Ha) Lpi^.^^dio Tj Sco ^'co FWHMco 
10^^ Lq K Jykms"^ 10^°K kms"^ pc^ kms~^ 



Mh 



M„ 



-dyn ^^^gas 

10^° Mq 10^ Mq 



16 36 55.04 


41 04 32.0 


2.189 2.186 2.192 


7.1±0.9 


>46 


0.34 ±0.07 


1.04±0.21 


13 12 36.01 


42 40 44.1 


2.240 2.243 


6.7±0.6 


>48 


<0.27" 


<0.71 


12 37 11.34 


62 32 31.0 


~1.996 ... 1.996 


16.7±0.7 


>45 


0.70 ±0.22 


1.93±0.60 



410±40 



>.4 sin^ 



8.3±1.7 

<5.6 
15.4±4.8 



a) For RG J131236, we set a limit for a 500 kms ^ FWHM line centered at the Hq redshift. 
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Fig. 3. — A comparison of the CO and radio luminosities for the 
three SFRGs described here, the SMGs from Greve et al. (2005) 
and Tacconi et al. (2006), the lensed LBGs detected in CO (Baker 
et al. 2004, Coppin et al. 2007, Kneib et al. 2005), the luminous 
2 ~ 2 BX/BzK galaxies undetected in CO (Tacconi et al. 2008) 
and the two z ~ 1.4 BzK galaxies detected in CO by Daddi et al. 
(2008). The solid line is the best-fitting relation with a form of log 
LJ-,Q = a log Lpipt+/3 to the local LIRGs and ULIRGs and the high- 
redshift SMGs from IGreve et al.l I I2005I1 (however it is no longer 
quite the best fit when radio luminosity is considered consistently 
across the populations). With gas conversions fixed, the SFRGs 
appear to have lower CO gas masses than SMGs, although if the 
radio-inferred SFRs are overestimated then the SFRGs could still 
lie on the plotted gas/SFR relation. 



of Lfir/Mh2 ~ 450 ± 17OL0Mq1 (IGreve et al.l 120051) . 
gas-to-dust mass ratios of ^ 200 (with a factor of a few 
uncertainty in the dust mass alone) and gas surface den - 
sities of cTgas -- SOOOM© yr'^ pc^ (|Tacconi et al.ll2006l ). 
Borys et al. (2005) estimate the average stellar mass for 
SMGs at z > 1.5 to be 3.2t?:6 x 10" M©, slightly larger 
than the average 1.8 x 10^^ M© for our three SFRGs. 

Overall, this comparison suggests that SFRGs may be 
somewhat smaller mass objects (lower stellar mass, lower 
CO mass and lower dynamical mass) than SMGs, but 
share with them a large radio luminosity. The SFRs 
of both SMGs and SFRGs are subject to sizable uncer- 
tainties, not least of which is the initial mass function 
(e.g., Baugh et al. 2005). Pope et al. (2006) have pointed 
out that SMG SFRs estimated from the 24 /im Spitzer- 
MIPS observations are lower than those estimated from 



the 850 /im or radio wavebands, although this could rep- 
resent an issue of relative calibrations of these indicators 
in this luminosity regime rather than intrinsic properties. 
We can also compare the SFRGs to local populations 
and less luminous star- forming galaxies at z ~ 2. Lo- 
cally, L'nn i ncreases with Lfir for (U)LIRGs, with the 
IGreve et al.l (|2005) sample of SMGs extending this trend 
out to the highest far-IR luminosities (> lO^'^L©). For 
comparison, in Fig. 3 we have plotted the SFRGs on the 
L^Q-LpiR, diagram along with SMGs lying on the local 
relation, as well as three LBGs from the literature within 
the considerable uncertainties in their far-IR luminosi- 
ties, three undetected BX/BzK galaxies (Tacconi et al. 
2008), and two CO-detected BzK galaxies (Daddi et al. 
2008 - these sources lie above the relation). The SFRGs 
lie somewhat below this relation. 

5. DISCUSSION 

As our calculations above have shown, it is very diffi- 
cult to estimate the precise gas masses for SFRGs due to 
various uncertainties. However, a potentially important 
result emerges from these CO observations: compared to 
SMGs, SFRGs appear to be significantly more efhcient 
at producing stars from a given molecular gas mass. If 
this is strictly true, then SFRGs cannot be interpreted 
as scaled up versions of local ULIRGs as Tacconi et al. 
(2006, 2008) have argued is the case for SMGs, since their 
gas masses appear to be lower than expected for their 
radio luminosities. There are two considerations to be 
taken into account here. Firstly, the far-IR luminosities 
and thus SFRs may be over-estimated from the radio, 
for instance if buried AGN were present. While none of 
these sources have AGN signatures in their UV or optical 
spectra, a deeply obscured AGN could still be driving a 
significant portion of the radio luminosity (e.g., Daddi et 
al. 2007b; Casey et al. 2008a). A further possible com- 
plication comes again from the radio selection of these 
objects. There is a ^--^0.25 dex scatter in the radio-FIR 
relation (Yun et al. 2001), and while locally there is no 
apparent correlation between SED shape and radio-FIR 
scaling, it is possible in our SFRGs that we are selecting 
galaxies which are amongst the lower 0.25 dex (weaker 
FIR per unit radio) . If the SFRs in our galaxies were sev- 
eral times lower then the efficiencies would be similar to 
the average SMG (Fig. 3). Secondly, the conversion from 
CO (3-2) to molecular gas mass may not be the same as 
for SMGs. If aco were greater than the '^l estimated for 
local ULIRGs (and as inferred to be correct for SMGs), 
these SFRGs could make up the shortfall in molecular gas 
mass from the average SMG (~ 4x), although the con- 
version would have to approach that typically adopted 
for the Milky Way aco = 4.6 (Solomon & van den Bout 
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Table 2 
Photometric Properties of SFRGs 



20cm 



850fim 
mJy 



350/im 



24/^m 



8.0/im 



5.8/im 



4.5^in 



3.6^in 

MJy 



/ R 

AB mag yl_B mag 



RGJ163655 43.9±7.1 -1.5±1.1 
RGJ131236 48.7±4.1 0.4±1.1 
RGJ123711 133.1±5.1 2.2±1.2 



13.7±6.9 < 150 14.4±2.7 

8.4±2.7 

531±8 37.5±0.8 



16.8±3.5 11.9±1.1 9.9±0.7 23.5 23.8 

< 11 7.0±1.1 5.4±0.7 24.2 25.0 

42.1±1.2 35.8±0.3 32.3±0.2 24.9 25.8 



Table 3 
Derived properties of the SFRGs 



SFE" 

L0/M0 



Mr^yr"^ Mr ^"■^'^ 



l0yi- 



SFR//Q 



SFR24fim ' 

Mpyr-i 



^gas 

Mq yr~^ kpc" 



M0 



gas /dust ' 



'^depletion 

(Myr) 



''gas to stars 

(Myr) 



< 255 600 70 1 X 10^^ 52 13.2 160 

< 1450 < 65 5 X 10"' - < 9.2 80 

880 ± 54 1330 35 4 x lO" 100 9.2 240 



RGJ163655 630 ± 120 46 ± 11 150 ± 50 
RGJ131236 610 ±91 33 ± 16 110 ± 40 
RGJ123711 1670 ±113 11 ± 10 16 ± 9 



a) SFR from 24fiin luminosity, assuming the average SED in Dale & Helou (2002). 

b) Star formation efficiency (SFE) is the SFR divided by the molecular gas mass. 

c) Stellar mass calculated as Borys et al. (2005), adopting their L^/Mq = 3.2. 

d) No limit is possible for RG J131236 since both gas and dust are limits. 



2005). This is unlikely given that the SFRGs often show 
clear evidence in high-resolution radio observations for 
merger-driven starbursts (Casey et al. 2008b). 

ft is noteworthy that while our observations have high- 
lighted an ultraluminous z ~ 2 population which may 
build stars in an extremely efficient mode (or at least as 
efficient as SMGs if their SFRs are overestimated by fac- 
tors of several), recent observations (Daddi et al. 2008) 
have identified a population of z ~ 1.5 galaxies which 
they detect in C0(2-l) exhibiting the opposite property: 
low-efficiency star formation. These Daddi et al. galax- 
ies are selected as large massive disk-like galaxies, and it 
is perhaps not surprising that they form stars in an ap- 
parently quiescent "spiral-galaxy" mode. Nonetheless, it 
is intriguing that galaxies in the high-redshift Universe 
have been discovered with such a wide range of star- 
forming efficiencies (from poor to extreme) all lying in 
the 10^^^^^-Lq regime. Massive galaxies are being built 
in a variety of modes in the z = 1 — 3 peak star- formation 
period. 

The short depletion times compared to the long times 
to form the stellar masses suggests we may be seeing 
SFRGs in the last phase of their current star-formation 
episodes. However, one would expect on average to find 
SFRGs as a population half-way through their gas con- 
sumption lifetimes. In this context, the small ages im- 
ply either a very high duty cycle or that the SFRs are 
over-estimated from the radio luminosity. We reiterate 
that the large FIR luminosity estimates for our SFRGs 
are based mainly on the radio/FIR relation which is ap- 
parently applies for SMGs (Kovacs et al. 2006), but is 
only marginally supported for SFRGs through the 24- /im 
luminosity and dust-corrected Ha measurements, ft is 
possible that we are over-estimating their star-formation 
activity, which would lead us to under-estimate the du- 
ration timescales above. 

We also note that the observed CO may be closely as- 
sociated with the star formation, and thus is probably 
warm and highly visible (i.e. with low aco)- It does 



not rule out the existence of a cooler less visible com- 
ponent (with high aco), not intimately associated with 
the current zone of star formation (e.g. in the inner disk), 
which may still become available within its own dynam- 
ical timescale to fuel star formation. The result above 
may therefore be affected by the selection for the highly 
visible component of the molecular gas. However, the 
cold gas would need to be fairly widely distributed (to 
ensure that it doesn't violate the dynamical limits on 
the total gas-|-stellar mass in the central regions), yet it 
must also be able to flow into the central regions on a 
timescale of <100 Myrs (near the limit of the gas sound 
speed ^ 100 km/s). 

In order to evaluate the implications of these results 
for massive galaxy formation, we recall that these three 
galaxies have typical characteristics for the SFRG pop- 
ulation at all wavelengths measured. Our results un- 
derline the major role of gas consumption over short 
timescales and with high efficiencies, characterizing rapid 
and strong merger-driven bursts as a major growth mode 
for both stellar mass and black holes in the distant Uni- 
verse. Even if the SFRs in these SFRGs were over- 
estimated by a factor of a few, they would remain 
ULIRG-class galaxies. If we assume that ~50% of the 
submm-faint SFRGs at z ^ 2 arc dominated by star for- 
mation at levels comparable to SMGs, we arrive at a 
density of ~ 5 x 10~^Mpc~'^, similar to that observed 
for SMGs (Chapman et al. 2003). Together, the SMGs 
and SFRGs represent a volume density 10 x smaller than 
measured for galaxies inferred to be forming stars at low 
efficiencies by Daddi et al. (2007b) which have space den- 
sities of order of lO'^Mpc^^. With SFRs a few to 10 x 
larger in the SMGs and SFRGs, the net effect is roughly 
equal numbers of stars being formed in both high- and 
low-efficiency modes at z ^ 2. 

Further study should ascertain whether other SFRGs 
follow a similar pattern to the galaxies studied in this 
paper. A substantial sample will allow the properties of 
the gas to dynamical mass ratio to be determined accu- 
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rately for the population, since we only currently have 
one well-detected line profile, and the dynamical mass 
determination is limited by the unknown inclination. 

6. CONCLUSIONS 

• We conclude that the radio luminosities of these 
SFRGs are higher for their overall mass (gas plus stel- 
lar) than for the SMGs (given that if the radio-SFRs are 
over-estimated for one class, they could well be for both). 

• We note that SMGs in general, and also these SFRGs, 
are outliers of the stellar mass-SFR correlation (Daddi 
et al. 2007a), probably due to the higher efficiency in 
forming stars for a similar stellar mass and CO luminos- 
ity. Lower gas masses in the SFRGs would imply even 
higher SFEs than the SMGs. By contrast, if the SFRs 
are significantly over-estimated by the radio, or the CO- 
to-H2 conversion were significantly different from SMGs, 
then SFRGs could have similar efficiencies to typical 
ULIRGs. Together with the apparent low-efficiency star- 
forming (U)LIRGs from Daddi et al. (2008), the SMGs 
and SFRGs with SFRs several to 10 x larger than the 
Daddi et al. galaxies suggest roughly equal numbers of 
stars being formed in both high- and low-efficiency modes 
at z ~ 2. Massive galaxies are being built in an impres- 
sive variety of modes in the z = 1 — 3 peak star- formation 
period. 

• If the radio-inferred SFRs are correct, then these 
SFRGs are more efficient star formers than SMGs, and 



cannot obviously be interpreted as scaled up versions of 
local ULIRGs as Tacconi et al. (2006, 2008) have argued 
is the case for SMGs. The SFRGs' radio luminosities 
are larger than would naturally scale from local ULIRGs 
given the gas masses or gas fractions. These observed 
gas masses and star-formation properties may be typical 
of the SFRG population and further work is justified to 
explore this population with improved statistics. 

• Our results underscore the fact that ultraluminous 
galaxies in the high-redshift Universe have been discov- 
ered with a wide range of star-forming efficiencies, the 
SFRGs apparently being one extreme. Massive galaxies 
are likely being built in a variety of modes in the z = 1 — 3 
peak star-formation period. 
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